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We study the localization properties of weakly interacting Bose gas in a quasiperiodic potential 
commonly known as Aubry-Andre model. Effect of interaction on localization is investigated by 
computing the ‘superfluid fraction’ and ‘inverse participation ratio’. For interacting Bosons the 
inverse participation ratio increases very slowly after the localization transition due to ‘multisite 
localization’ of the wave function. We also study the localization in Aubry-Andre model using an 
alternative approach of classical dynamical map, where the localization is manifested by chaotic clas¬ 
sical dynamics. For weakly interacting Bose gas, Bogoliubov quasiparticle spectrum and condensate 
fraction are calculated in order to study the loss of coherence with increasing disorder strength. 
Finally we discuss the effect of trapping potential on localization of matter wave. 
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I. INTRODUCTION 

In recent years ‘Anderson localization’ of particles and 
waves has regained interest in quantum many particle 
physics. ‘Anderson localization’ is a remarkable quantum 
phenomenon for which the propagating wave becomes 
exponentially localized in the presence of disorder (TJ |2] 
and is a well studied subject in the context of electronic 
systems in presence of disorder. Like particles, waves 
can also localize in disordered medium. Localization 
of ‘matter wave’ has recently been observed in exper¬ 
iments on ultracold Bose gases in presence of speckle 
potential [3] and in bichromatic optical lattices [4]. For 
certain parameters, Bosons in bichromatic lattice can be 
mapped on to Aubry-Andre (AA) model [5] with ‘quasi¬ 
periodic’ potential. Also localization of light in AA 
model has been observed experimentally [6]. Unlike An¬ 
derson model with random disorder, the AA model ex¬ 
hibits localization transition in one dimension at a criti¬ 
cal strength of the potential. Apart from that the quasi¬ 
periodicity of AA model gives rise to various interest¬ 
ing spectral properties [7]. In recent experiment [8] ef¬ 
fect of interaction on AA-model has been studied. Re¬ 
pulsive interaction plays an important role in the for¬ 
mation of correlated phases like ‘Bose-glass’ phase 
[12] . Also the dynamics and diffusion of interacting Bose 
gas in presence of disorder have been investigated both 
experimentally [131 Ej nnd theoretically [T^ [TB] , In re¬ 
cent years ‘many body localization’ m and localization 
at finite temperature [T8ll2Q] have generated an impetus 
to study disordered Bose gas. 

In this work we investigate localization of weakly in¬ 
teracting Bose gas in AA potential. Apart from calcu¬ 
lating the ground state properties, we also study local¬ 
ization of wavefunction by using a classical dynamical 
map. The paper is organized as follows: In section II, 
we discuss non-interacting Bose gas in quasiperiodic po¬ 
tential. In subsection A, we review the AA model and 
discuss the single particle localization properties. Local¬ 
ization in non-interacting system using dynamical map 
approach is presented in subsection B. In section III, 


various physical quantities like inverse participation ra¬ 
tio and superfluid fraction are computed to investigate 
localization of weakly interacting Bose gas within mean- 
field theory. Effect of non-linearity due to interactions 
on the dynamical map is studied. Further, we compute 
the Bogoliubov quasiparticle energies and amplitudes to 
investigate quantum fluctuations. We study the localiza¬ 
tion of both Bogoliubov amplitudes and non-condensate 
densities. 

Effect of trapping potential on localization and effect 
of disorder on the center of mass motion are presented in 
section IV. Finally we summarize our results in section 

V. 


II. NON-INTERACTING BOSONS IN 
QUASIPERIODIC POTENTIAL 

In the original experiment [4] localization of ultracold 
Bosons has been studied using bichromatic optical lattice 
which can be mapped on to Aubry-Andre model within 
tight binding approximation and for certain parameter 
regime [21]. The Aubry-Andre model is defined by the 
Hamiltonian, 

H = -J y~^(|n)(n+l| + |n+l)(n|)+A E cos(27r/3n)|n)(n| 

n n 

( 1 ) 

where |n) is a Wannier state at lattice site n, J is nearest 
neighbour hopping strength, A is the strength of onsite 
potential and period of potential is determined by /d. In 
the rest of the paper, we would be working in the units 
in which J = 1. 

For A = 2, the Hamiltonian 0 is equivalent to the 
Harper model [22] describing the motion of an electron 
in a square lattice in the presence of a perpendicular 
magnetic field, where the flux 4> through each plaque- 
tte in units of flux quantum 4>o = hje is given by the 
parameter (3 = 4>/4>o. Energy spectrum of this problem 
gives rise to the well known ‘Hofstadter butterfly’ [23]. 
The Hamiltonian given in Eq.Q, poses very interesting 
properties for irrational values of /3. When (3 is chosen 
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to be a ‘Diophantine number’, A A model undergoes a 
localization transition at a critical value of the potential 
strength A = 2|5l[2l]. On the contrary, localization tran¬ 
sition is absent in one dimensional Anderson model with 
random disorder. The quasiperiodic potential generates 
a correlated disorder in the AA model. 


A. Single particle localization 


In order to study localization transition in AA model, 
we choose P = {^/b — l)/2 which is the inverse of ‘golden 
mean’ and a Diophantine number . This is particularly 
helpful since rational approximation of P can be done 
by Fibonacci series. Although incommensurability of the 
potential with the underlying lattice plays an crucial role 
in AA model, for numerical studies one can approximate 
P = Fn-ijFn^ where F^ is nth Fibonacci number for 
sufficiently large value n. This rational approximation 
of P fixes the lattice size Ng = Fn to impose periodic 
boundary condition [25]. 

Duality in AA model can be shown by introducing new 
basis states in momentum space. 


\k) Ng E exp(i27r/C/dn)|n). (2) 

n 

The dual model is obtained by substituting Eq@ in Eq0 



\ k){k F l \ h.c + ^ cos{27rPk)\k){k\ 
. k k 


( 3 ) 

It is important to note that AA model becomes self-dual 
at a critical coupling A = 2, where the localization transi¬ 
tion occurs. To obtain the eigenvalues and eigenfunctions 
of the Hamiltonian we expand the state vector in terms of 
Wannier states, where is the wave- 

function at nth lattice site. The eigenvalue equation of 
the Hamiltonian Eqj^is reduced to a discrete Schrodinger 
equation. 


- (V’n+I + V’n-i) + Acos(27r/3n)'*/)„ = (4) 

where e is the energy eigenvalue. The degree of local¬ 
ization of a normalized state Ipj) can be quantified by 
‘inverse participation ratio’ (IPR) /, 

i = ( 5 ) 

n 

The wavefunction of an extremely localized particle at a 
site no is given by ipn = ^n,no , for which the IPR becomes 
unity. On the other hand, for the completely delocalized 
wavefunction IPR is l/Ng which vanishes 

in the thermodynamic limit. For A A model all energy 
eigenfunctions in real space are exponentially localized 
and IPR sharply increases to unity above the critical 
coupling A = 2. Due to the duality of A A model, the 
localization of wavefunction in real space and in dual mo¬ 
mentum space shows opposite behavior. In real space the 
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FIG. 1: IPR of ground state wavefunction as a function of 
strength of the potential A. IPR in real space F and in dual 
momentum space A are shown by solid (black) line and dashed 
(blue) line respectively. 


wavefunctions are localized above A = 2, wheras localiza¬ 
tion in dual momentum space occurs for A < 2. The IPR 
of the ground state wavefunction in real space and in dual 
monentum space as a function of A are shown in Fig{^ 
It is interesting to note that IPR in real and momentum 
space intersects at the self dual point A = 2. 

The spectrum of AA model also shows various interest¬ 
ing properties. For A = 0 it has usual single band energy 
spectrum of periodic lattice. Addition of the potential 
generates quasiperiodic structure and destroys the band 
like dispersion. Successive rational approximation of P 
generates more periodicity over the underlying lattice, 
which in turn breaks the original energy band into many 
subbands and it leads to opening of band gaps. The 
energy spectrum of AA model is obtained by numerical 
diagonalization and is depicted in Figj^ for increasing 
values of A. The variation of energy gaps with the cou¬ 
pling strength can be noticed in this figure. The spectral 
statistics and distribution of energy gaps are analyzed in 
[ 7 ] . Mathematically it can be shown that the energy spec¬ 
trum of AA model forms a Cantor set [26]. Self-similarity 
in energy levels can be understood from the integrated 
level density. 


7V(e) = ^0(e-e,), (6) 

i 


where is the i’th eigenvalue and 0{x) is heaviside step 
function. The integrated density of states N{e) shows 
‘Devil’s staircase’ like fractal structure, which is evident 
from Figj^ where we plotted the normalized integrated 
density of states N{e) as a function of scaled energy 
within the interval of zero to one. For different values 
of A, the Devil’s staircase structures of N{e) do not over¬ 
lap, which indicates that the fractal dimension changes 
with the strength of the potential A. 

Transport properties also show significant changes in 
localization transition. Due to spatial localization of sin¬ 
gle particle states transport coefficients vanish in the 
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FIG. 2: (a) Eigenenergies of A A model as a function of A 
for system size Ng = 144. (b) Devil’s staircase structure of 
the integrated density of states N{e) scaled by N{emax) as a 
function of scaled energy = (e - emin)l{emax - Cmm) at 
A = 2 and for different system sizes. Here Cmin and Cmax are 
the minimum and maximum energy eigenvalues. 


thermodynamic limit. In electronic systems the conduc¬ 
tivity vanishes exponentially with the system size due 
to ‘Anderson localization’. For neutral superfluid corre¬ 
sponding physical quantity is ‘superfluid fraction’(SFF), 
which is measured by generating a superflow by apply¬ 
ing a phase twist [27]. In presence of the phase twist the 
original Hamiltonian in Eqa becomes, 

Ns 

Hq = — |n)(n + l| -|-h.c) + A cos(27r/3n)|n)(n|, 

n=l 

(7) 

where 0 is arbitrarily small phase difference across the 
boundary. For one dimensional system with periodic 
boundary condition this is exactly similar to a quantum 
ring in presence of a flux which generates supercurrent 
through the ring. The superfluid fraction fs is defined 
as [28]. 


( 8 ) 

where Eo{&) is the ground state energy of the Hamilto¬ 
nian Hq with arbitrary small value of 0, and F^o(O) is the 
ground state energy of the original Hamiltonian given in 
EqUl To obtain the ground state energy upto 0^ order, 
the Hamiltonian Hq is expanded as. 


Hq 


n/ 2Nr 


(9) 


where we define a current operator J = + 

1| — h.c) and the usual kinetic energy T = — X]n(l^)(^ + 
1| + h.c). Using second order perturbation, the SFF in 
Eqnj^ can be written as. 


fs — ~^('0oi^i'0o)+ 


~ ^0 


( 10 ) 


where e^, \ 2 pi) are eigenvalues and eigenfunctions of 
Hamiltonian given in Eqj^and |'0o) is ground state wave- 
function. Variation of SFF with the coupling strength of 





FIG. 3: Superfluid fraction fs (solid line) and IPR of ground 
state in real space Ir (dashed line) as a function of A. 


the potential is shown in Figj^ For increasing values of A 
the SFF decreases from unity and vanishes at the critical 
point A = 2. 


B. Hamiltonian map 

Localization phenomena can also be studied by an al¬ 
ternative method of classical Hamiltonian map(CHM), 
which is very useful to obtain analytical estimate of lo¬ 
calization length [2^ l3Q] . The discrete Schrodinger equa¬ 
tion given in Eqj^can be written as following dynamical 
map. 


Pn+I = Pn + (Acos(27r/3n) - e - 2)Xn, (11) 

^n+l ~ T Pn+l (1^) 


where the classical dynamical variables are given by, Xn = 
and Pn = 'f’n — Here the wavefunction 'ipn 

plays the role of position in CHM and the number of 
lattice site becomes the number of iteration or time axis 
of the dynamics. We can choose initial real wavefunctions 
-00 and 'ijji (or equivalently Xi and pi) and evolve the 
dynamical variables by the transfer matrix. 


(Acos(27r/3n) — e — 1) 1 
(A cos(27r;5n) — e — 2) 1 


(13) 


By iteration of this map we obtain the asymptotic behav¬ 
ior of the dynamical variables for a fixed energy eigen¬ 
value e. Mathematically it can be shown that the dy¬ 
namical instability of the above Hamiltonian map cor¬ 
responds to the localization phenomena. Exponentially 
localized wavefunctions asymptotically fall off as ^ 
where ^ is the localization length. But for ar¬ 
bitrary choice of the initial wavefunction, the solution 
of Eq0 also has an exponentially growing component 
^ ^ Hence the localization is manifested by the expo¬ 

nential growth of the dynamical variables of the Hamilto¬ 
nian map which gives rise to chaos in the classical phase 
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FIG. 4: Phase portraits of CHM corresponding to Eq[T^with 
a fixed energy eigenvalue close to the middle of the energy 
band and for an ensemble of initial conditions. Fig. (a) to (d) 
corresponds to A=0, 0.6, 1.2, 1.8 respectively. 


space. Whereas the extended states represents periodic 
motion in the phase space. To understand the localiza¬ 
tion transition in AA model, we choose an eigenvalue 
close to the band center and calculate the phase space 
trajectories of the Hamiltonian map for an ensemble of 
random initial values of x and p. The phase portraits of 
EqfT^ are shown in Figj^ for increasing values of A. In 
the delocalized regime with small disorder strength A, the 
phase portrait consists of closed elliptic orbitals as seen 
from Fig.3a. Increasing the strength of quasiperiodic po¬ 
tential A leads to diffusive behavior at the outer part of 
the phase space and regular portion of phase space with 
periodic orbits is reduced (see Fig.3b and Fig.3c). Fi¬ 
nally for A close to the critical value, an instability sets 
in the CHM and phase space trajectories in Fig.3d show 
chaotic behavior which indicates localization transition. 

In dynamical systems ‘Lyapunov exponent’ (LE) is a 
measure to quantify chaos. Since the number of lattice 
site plays the role of time in the classical map, the LE 
corresponds to the inverse localization length of the wave- 
function. For periodic motion LE vanishes and in the 
chaotic regime the non-vanishing LE gives finite localiza¬ 
tion length To calculate the LE we construct the ma¬ 
trix Uns = Tn^Tn ^-1 . T 2 T 1 by multiplying the transfer 

matrices sequentially for Ng iteration. The LE I can 
be obtained by using the formula. 


I = 


lim 

Ng —^00 


log(A+) 

Ns 


(14) 


where A+ is the largest eigenvalue of the matrix U 
LE is obtained from ^ = 1/1. The localization length can 


also be calculated by using Thouless formula [31]. Us¬ 
ing the duality of AA model the analytical expression of 
LE is given by ^ = log(A/2)[5|, which is independent of 
energy. Within CHM approach we numerically compute 
the localization length ^ from Eq(g using QR decom¬ 
position method [32]. Numerically obtained localization 
length as a function of A for different energy eigenvalues 
are compared with the analytical result in Figj^ 
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FIG. 5: Variation of Lyapunov exponent or inverse localiza¬ 
tion length l/(f with disorder strength A. Dotted and dot- 
dashed lines represent numerically obtained Lyapunov expo¬ 
nents from Eq[g and for comparison the analytical result is 
shown by solid (red) line. 


For the localization transition in AA model with 
quasiperiodic potential, the choice of /3 as irrational num¬ 
ber (particularly a Diophantine number) plays a crucial 
role m- To understand this mathematical condition, we 
analyze the phase space of CHM for successive rational 
approximation of (3^ which is shown in Figj^ Using Fi¬ 
bonacci series, the rational approximation of {3 can be 
written as /3n = Fn-i/Fn and the potential has period¬ 
icity Fn for successive integer n. For sufficiently large 
value of n, [3n approaches to the inverse of ‘golden mean’ 
and the potential becomes quasiperiodic. In the local¬ 
ized regime, for A = 2.2 the phase space trajectories of 
CHM with increasing order of rational approximation of 
P are presented in FigJ^ (a) to (c) for fixed initial condi¬ 
tion. Even in the localized regime phase space contains 
periodic orbits for As shown in Fig.6b to Fig.6d, the 
periodic orbits break as p approaches to the inverse of 
‘golden mean’ by successive rational approximation and 
finally dynamics become chaotic which is consistent with 
the localization phenomena. 
























5 



(c) 


(d) 


FIG. 6: Phase space trajectories with successive rational ap¬ 
proximation of /3 and for A = 2.2. Here (a) to (c) corresponds 
to = 0.6, ydr = 0.615, ydi 2 = 0.618 respectively with same 
initial conditions and (d) corresponds to /3i4 = 0.618 with an 
ensemble of initial conditions. 


III. LOCALIZATION OF WEAKLY 
INTERACTING BOSE GAS 

Interacting Bosons in presence of qnasiperiodic poten¬ 
tial can be described by Bose-Hnbbard model [9], 

H = — y^(atajH-h.c)+A cos(27r/3i)ni+-^ ni{ni — l) 

(ij) i i 

(15) 

where at(ai) are creation (annihilation) operators for 
Bosons at site hi = ata^, and U is the strength 
of onsite repulsive interaction. Above quantum many- 
body Hamiltonian can capture various correlated phases 
of strongly interacting Bosons which undergoes quan¬ 
tum phase transition. It is known that Bose-Hubbard 
model in the presence of random disorder can give rise 
to ‘Bose glass’ phase [9]. For sufficiently weak interac¬ 
tion strength U and for high average density of Bosons a 
‘quasi-condensate’ may form in one dimensional systems 
[33] . In this regime, one can replace the bosonic operators 
ai by a classical field which represents the macroscopic 
wavefunction of the ‘quasi-condensate’. Minimization of 
the classical energy corresponding to Eqm leads to the 
‘discrete nonlinear Schrodinger equation’ (DNLS), 

-{'ipi+i+ipi-i)+Xcos{2'Kl3i)xl)i + U\xl)i\'^ipi = iiipi, (16) 

where /i is the chemical potential, and the normaliza¬ 
tion of the wavefunction of number of Bosons gives 

We numerically find out the ground state wavefunc¬ 
tion of Eq{l^ by imaginary time propagation method. 



(a) 



FIG. 7: (a)Inverse Participation Ratio(I) of condensate wave- 
function with increasing A for different UNb. (b) Spatial vari¬ 
ation of condensate density with A for UNb = 0.05, where i 
denotes lattice site index. 


For weak disorder the wavefunction is extended and nu¬ 
merical convergence is fast, whereas close to the localiza¬ 
tion transition more care is needed to obtain the ground 
state since many metastable states appear in this region. 
To obtain the degree of localization of the ground state 
wavefunction in presence of repulsion we calculate the 
IPR in real space using Eqj^ The variation of IPR of 
condensate wavefunction with increasing strength of dis¬ 
order A for different values of effective interaction U is 
shown in FigJ^a). The IPR becomes nonzero for A > 2 
and increases with a slower rate compared to the non¬ 
interacting system, indicating spreading of the wavefunc¬ 
tion due to the repulsive interaction. The degree of lo¬ 
calization decreases with increasing strength of repulsive 
interaction UNb. In Fig.7(b) the spatial variation of the 
wavefunction with disorder strength A is represented by 
color scale plot. It is evident from this figure that sin¬ 
gle site localization is not favorable energetically and the 
wavefunction is localized at almost degenerate but spa- 
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tially separated sites. This particular feature of the frag¬ 
mented condensate has also been studied in [34]. Due 
to the multisite localization the IPR is much less than 
unity and increases very slowly with A. A change in the 
slope of IPR with increasing A occurs when the number 
of localized sites decreases and we have checked that fi¬ 
nally at very large value of A the wavefunction becomes 
localized at a single site. 

To study the interplay between disorder and interac¬ 
tion in the transport properties of dilute Bose gas, we cal¬ 
culate the superfluid fraction fs . To generate a superflow 
we introduce a small amount of phase twist (0 ^ 0.1) in 
the hopping term of DNLS (in Eqp^ similar to Eqj^ 
and then the superfluid fraction is computed using the 
formula, 


. _ ,,2^d(0)-i^d(O) 

•'« ~ /a2 


(17) 


where Eci{Q) is the classical energy corresponding to the 
Hamiltonian in EqfTs] 


Eci{e) 


“E +h.c 

{ij) 



A cos(27r/3i)|0i|^ 

i 


(18) 


ten in the form of nonlinear classical map, 

Pi+i = Pi + {\ cos{2tt j3i) - ii-2)xi + UNbx\,{lQ) 

Xi+i = Xi+ pi+i ( 20 ) 


where, Xi = c^i and Pi = (pi — <pi-i- The repulsive 




(c) 


(d) 


The wavefunction (pi minimize Eci{&) and is normalized 
to unity. The superfluid fraction fs as a function of A for 
different values of repulsive interactions UNij is shown in 
Eigj^ The SEE of weakly interacting Bose gas obtained 
from ‘density matrix renormalization group’ also shows 
similar behavior [35|. Due to the repulsive interaction, 
fs vanishes at larger strength of quasiperiodic potential 
A > 2, however the IPR rises from zero at A ~ 2. This 
behavior is different from the noninteracting AA model. 



FIG. 8: Variation of ‘superfluid fraction’ fs with A for differ¬ 
ent interaction strength UNb. 


FIG. 9: Phase portraits of GHM for increasing values of A 
and with interaction strength UNb — 0.5. Figure (a) to (d) 
correspond to A=0, 0.2, 0.9, 1.8 respectively. 

nonlinearity gives rise to an unstable classical potential 
^ due to which the classical trajectories be¬ 

come unstable at large values of phase space variables. 
To avoid this problem we choose relatively small region 
of phase space within which the potential is metastable, 
and study the effect of disorder in the phase space dy¬ 
namics. In Figj^ for a fixed value of nonlinearity UN^ 
we show the phase portrait with an ensemble of initial 
configurations for increasing values of disorder strength 
A. Similar qualitative features like non interacting sys¬ 
tem are also seen in the phase space dynamics of DNLS, 
however the classical periodic orbits are modified due to 
the nonlinearity. 

A. Fluctuations within Bogoliubov approximation 

So far we studied the weakly interacting Bosons within 
mean-field approximation using macroscopic wavefunc¬ 
tion for the condensate. It is also important to analyze 
the quantum fluctuations induced by the quasiperiodic 
disorder potential. Within Bogoliubov approximation, 
the quantum field operators can be approximated by. 


We also investigated the localization properties of 
DNLS by Hamiltonian map approach. Eqj^can be writ- 










, ( 21 ) 
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where ijji is the macroscopic wavefunction of the conden¬ 
sate satisfies Eqj^ are amplitudes corresponding 

to z/th eigenmode with bosonic operators Bogoli- 

ubov quasiparticle energies can be obtained from, 

-{Ui+i + Ui-i) + [Acos(27r/3i) + 2U\'4>i\^ - ii]ui 

+U'tpiVi = coui, ( 22 ) 

-{vi+i +Vi_i) + [Acos(27r^i) + 2J7|V’iP - 

= -wvj, (23) 

and normalization condition gives ) = 

Suu'‘ 





140 

110 

80 

50 

20 


0 12 3 4 

A 

(c) 

FIG. 11: Variation of (a) condensate fraction Nc/Nb and (b) 
inverse participation ratio / of non-condensate density pnc 
with A for different interaction strengths UNb as shown in the 
figure, (c) Spatial distribution of pnc with increasing disorder 
strength A keeping UNb = 1. For all the above figures, Ng = 
144 and Nb = 500. 



FIG. 10: (a) Bogoluibov excitation spectra with increasing A 
for UNb = 1. (b) Devil’s staircase structure of the integrated 
density of states N{uj) (scaled by N{(jUniax)) as a function of 
scaled energy Ug = {uj - ojmin)/(cjrnax - (^min) for A = 2.5 
and UNb = 1 . Here Uniin and Umax are minimum and max¬ 
imum energy eigenvalues, respectively, (c)-(d) : Variation of 
normalized density of Bogoliubov amplitudes and 
respectively with A. The index ly represents an energy state 
at the center of the band with Uu ^ 2. 

For a given interaction strength U and average density 
of Bosons no = N^/Ns, we first numerically obtain the 
ground state macroscopic wavefunction then diago¬ 
nalize Eqj^to calculate Bogoliubov quasiparticle ener¬ 
gies and amplitudes . The Bogoliubov energy spec¬ 

trum for increasing strength of disorder A is depicted in 
Fig@ a) for a fixed value of interaction UN^ = 1 and 
Ns = 144. The normalized integrated density of states 
N{ijj) also shows Devil’s staircase like structure as shown 
in Figp^b). Due to the quantum fluctuations, depletion 
of the condensate occurs and the non-condensate den¬ 
sity pnc at zero temperature can be obtained from the 
Bogoliubov theory, 

/0nc{*) = EKl'- (24) 


The condensate fraction is given by A^c/^6 = 1 “ 
'Thiv dimensional system the noncon¬ 

densate fraction diverges as log(A^s), which prohibits the 
formation of condensate in the thermodynamic limit. 
However for sufficiently weak interaction a quasi con¬ 
densate can form in quasi one dimensional and finite 
system [33]. To investigate the disorder induced quan¬ 
tum fluctuation in a quasi-condensate in a finite lattice 
with Ns = 144, we calculate the condensate fraction 
with increasing disorder strength A, which is shown in 
Figl^ a). In absence of disorder the quantum depletion 
is small in weakly interacting gas of Bosons and increases 
with the interaction strength. With the increase in dis¬ 
order strength A the condensate fraction remains close to 
unity in the delocalized regime and then rapidly decreases 
around the critical point A ~ 2. This qualitative feature 
(as shown in FigjTTj^a)) indicates that enhanced quantum 
fluctuations near the localization transition can destroy 
the quasi condensate above A ~ 2 and strongly correlated 
phases can appear. The variation of non-condensate den¬ 
sity Pnc with disorder strength A is depicted in FigjTTj^c) 
by color scale plot. It is interesting to note that the 
distinct feature of multisite localization for A > 2 is ob¬ 
served even for non-condensate density. Also the IPR of 

















































normalized non-condensate density shows similar behav¬ 
ior as that of the condensate wavefunction and increases 
from A ~ 2 (as seen from Fig{^ c)). Localization of 
Bogoliubov quasiparticles and enhancement of quantum 
fluctuations in presence of quasiperiodic potential partic¬ 
ularly for A > 2 are clearly evident from this analysis. 


IV. LOCALIZATION IN THE PRESENCE OF A 
TRAP 

Although the localization transition in AA model oc¬ 
curs in thermodynamic limit, in real experimental setup a 
weak trapping potential is always present in order to con¬ 
fine the ultracold atoms. Main features of the localiza¬ 
tion can also be observed in trapped system provided the 
length scale of the trapping potential is larger compared 
to the localization length. Additionally some interesting 
effects due to the trap can also be seen. A harmonic trap 
is introduced by adding a potential Vi = ~ 

in the Hamiltonian (Eqjl^, where ujho is related to the 
trapping frequency, i is the site index, and ic is the cen¬ 
ter of the trap. Although in a trap the wavefunction is 



A 

(a) 



(b) 

FIG. 12: (a) Variation of inverse participation ratio / with 
increasing A for different trap frequencies uho (as shown in 
figure) keeping UNh = 0.5 and Nh = 500. (b) Density dis¬ 
tribution of harmonically trapped condensate for different A 
and for UNb = S, Nb = 500 and ojho = 0.05. 

always localized, but in a weak trapping potential the 
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FIG. 13: (a)Genter of mass motion for different disorder 
strength A and for ujho = 0.1, UNb = 2.5 and Nb = 500. For 
higher A oscillation becomes incoherent.(b)-Goherence factor 
which is measure of overall coherence in the system decreases 
for increasing A. 


width of the wavepacket is sharply reduced due to dis¬ 
order induced localization. The density profile of the 
condensate in presence of a harmonic trap for different 
disorder strength A is shown in Fig(T^. We calculate 
the IPR of the ground state wavefunction with increas¬ 
ing disorder strength A and an increase of IPR around 
A ~ 2 is observed as expected. However the IPR does 
not vanish in the delocalized regime A < 2 and takes a 
small value due to the tapping potential. We have noticed 
earlier that in the absence of a trap the IPR increases 
very slowly after the localization transition due to the 
repulsive interaction and the wavefunction is localized at 
spatially separated sites with quasi-degenerate energies. 
These quasi-degeneracy of onsite energies can be lifted 
by introducing a harmonic trap which leads to enhance¬ 
ment of the degree of localization which is elucidated in 
FigfT^ where a rapid increase of IPR to unity is shown 
by increasing the trap frequency by a small amount. 

Collective oscillation of a Bose-Einstein condensate in 
a trap can also show its superfluid properties. A small 
displacement of the condensate from the center of the 
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trap generates center of mass (COM) oscillation which is 
a well studied collective mode of Bose-Einstein conden¬ 
sate. Here we numerically study the motion of COM of 
trapped condensate in presence of quasi periodic disor¬ 
der. Since the superfluidity is affected by the disorder, 
we calculate the coherence factor ^ of the oscillating con¬ 
densate which is given by [36] 

^ = (25) 

i 

We can notice that the above expression contains infor¬ 
mation of relative phase difference between neighboring 
sites and gives a quantitative measure of overall 
phase coherence of the oscillating condensate. 

In Figp^ (a) and (b), we have shown the COM motion 
of the trapped condensate and coherence factor of corre¬ 
sponding macroscopic wavefunction for increasing disor¬ 
der strength A. It is clear that the coherence of the time 
dependent wavefunction decreases with increasing disor¬ 
der. It is also interesting to study the variation of conden¬ 
sate fraction with disorder. In one dimensional trapped 
condensate the divergence of non-condensate fraction is 
less severe due to the finiteness of the trap. Like the 
homogeneous system, we expect enhancement of quan¬ 
tum fluctuation due to disorder which may destroy the 
condensate around a critical disorder strength A 2. 

V. CONCLUSIONS 

In conclusion, we have investigated localization of 
both non-interacting Bosons as well as weakly interacting 
quasi-condensate in presence of a quasiperiodic potential. 
Apart from calculating various physical properties, un¬ 
derstanding localization transition through the approach 
of classical Hamiltonian map is one of the main result of 
this work. 

In the non-interacting Aubry-Andre model the local¬ 
ization transition can be identified by the vanishing ‘su¬ 
perfluid fraction’ and the rise of IPR at the critical 
strength of quasiperiodic potential A = 2. A classical 
Hamiltonian map is constructed from the Schrodinger 
equation. The phase space trajectories of the correspond¬ 
ing classical map shows periodic orbits for small disor¬ 
der strength. With the increasing disorder strength A 
chaotic behavior is observed at the outer region of the 
phase space and finally near the critical value A = 2 all 
periodic orbits are destroyed due to the onset of chaos 


in the localized regime. The parameter [3 = {^/b — l)/2 
is chosen to be an irrational Diophantine number (inverse 
of golden mean) which is an essential requirement for lo¬ 
calization transition. This has been elucidated by means 
of Hamiltonian map for A > 2 where the destruction of 
periodic orbits by successive rational approximation of /3 
indicates onset of localization. 

We have investigated the localization of weakly inter¬ 
acting Bose gas in quasiperiodic potential by mean-field 
approach. Unlike non-interacting case, the vanishing of 
SEE and rise of IPR does not take place at same strength 
of disorder A. Due to the repulsive interactions, the 
macroscopic wavefunction is localized at many sites with 
quasi-degenerate energies for A > 2. The multisite local¬ 
ization of the interacting system is manifested by a much 
slower increase of IPR starting from A = 2. With increas¬ 
ing the strength of quasiperiodic potential, the number of 
sites over which the condensate is localized decreases and 
finally the wavefunction becomes localized at a single site 
for very large value of A and IPR approaches to unity. 
The SEE decreases with increasing disorder strength A 
and vanishes at A > 2 due to the repulsive interaction. 
The repulsive interaction gives rise to an unstable non¬ 
linear potential in the CHM approach, due to which the 
stable region of phase space decreases. In the phase por¬ 
trait the stable region containing the periodic orbits de¬ 
creases with increasing A and finally the onset of chaos at 
A ~ 2 signifies localization of the wavefunction. Further, 
the Bogoliubov quasiparticle spectrum has been calcu¬ 
lated numerically. We notice that disorder enhances the 
quantum fluctuations due to which the condensate frac¬ 
tion of a finite system decreases rapidly around A ~ 2. 
This indicates the possible formation of glassy phase and 
multisite localized insulators. 

Finally we also considered the effect of trapping poten¬ 
tial on the localization transition. In the localized regime, 
the number of sites over which the wavefunction is local¬ 
ized reduces due to the presence of a trap which has been 
shown from the rapid increase in the IPR by tuning the 
trap frequency. The center of mass motion of the con¬ 
densate in a harmonic trap also shows the signature of 
localization. The center of mass oscillations become in¬ 
coherent with increasing disorder. To summarize, the 
present study provides a clear picture of localization of 
non-interacting and weakly interacting Bose gas in the 
presence of a quasiperiodic disorder and it reveals var¬ 
ious interesting features which are interesting for both 
academic point of view, as well for future experiments. 
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